arXiv: 1503.0683lv2 [astro-ph.CO] 20 Oct 2015 


Limits of Quasi-Static Approximation in Modified-Gravity Cosmologies 


Ignacy Sawicki 

Department de physique theorique, Universite de Geneve, 
Quai Ernest Ansermet 24, 1211 Geneve, Switzerland* 


Emilio Bellini 

Institut de Ciences del Cosmos, 

Universitat de Barcelona, IEEC-UB, 

Marti i Franque 1, E08028 Barcelona, Spain * 

We investigate the limits of applicability of the quasi-static approximation in cosmologies fea¬ 
turing general models of dark energy or modified gravity. We show that, at best, the quasi-static 
approximation breaks down outside of the sound horizon of the dark-energy, rather than the cos¬ 
mological horizon as is frequently assumed. When the sound speed of dark energy is significantly 
below that of light, the quasi-static limit is only valid in a limited range of observable scales and 
this must be taken into account when computing effects on observations in such models. As an 
order of magnitude estimate, in the analysis of data from today’s weak-lensing and peculiar-velocity 
surveys, dark energy can be modelled as quasi-static only if the sound speed is larger than order 1% 
of that of light. In upcoming surveys, such as Euclid, it should only be used when the sound speed 
exceeds around 10% of the speed of light. In the analysis of the cosmic microwave background, 
the quasi-static limit should never be used for the integrated Sachs-Wolf effect and for lensing only 
when the sound speed exceeds 10% of the speed of light. 


I. INTRODUCTION 

Dynamical dark energy (DE) and modified gravity 
(MG), modelling the observed acceleration of late-time 
cosmological expansion, usually require a new degree of 
freedom beyond those present in standard A-cold dark 
matter (ACDM) cosmology. 1 The full dynamics of this 
new degree of freedom in general models of DE/MG are 
usually quite complicated. However, since galaxy-survey 
data are still mostly available only on scales small com¬ 
pared to the cosmological horizon, the quasi-static (QS) 
approximation is frequently used to approximate the full 
DE/MG dynamics and interpret observations. Roughly, 
this amounts to neglecting terms involving time derivat¬ 
ives in the Einstein equations for perturbations and only 
keeping spatial derivatives. Such a procedure is thought 
to be valid on sufficiently small scales, typically assumed 
to be those well inside the cosmological horizon. It is well 
known from numerical studies that in the case of quint¬ 
essence [4, 5], f(R) gravity [6] or the covariant galileon 
models [7, 8] this approximation is good enough on ob¬ 
servable linear scales (see e.g. [9-11] and in particular ref. 
[12]). On the other hand, general perfect-fluid DE can 
show some discrepancy between the full and QS solutions 
at late times [13]. On non-linear scales, N-body simula¬ 
tions have been used to study symmetron-screened [14] 
non-minimally coupled quintessence [15, 16] and f(R) 
gravity [17] with largely the same conclusion: the effects 
of the time-derivative terms are small. But how general 


1 Although modifications of general-relativistic (GR) constraints 
also exist, e.g. refs [1—3]. 


are these results? 

The solutions in the dust-dark-energy dynamical sys¬ 
tem which models the late universe can indeed always be 
described using two functions of space and time [18-25] 
without reference to the extra degree of freedom. Tak¬ 
ing the QS limit allows one to remove the dependence 
on the initial conditions for the DE/MG and to reinter¬ 
pret what is a parameterisation of a particular solution 
as a universal description of a model. In the extreme 
QS limit, all the scale-dependence is also neglected and 
the DE model is described using just two functions of 
time: the effective Newton’s constant p, and the grav¬ 
itational slip parameter. Allowing for scale dependence 
in principle gives a closer description on a wider set of 
scales [26-28], but the question arises to what extent all 
of these approximations are valid. 

Recently, effective-field-theory (EFT) methods [29-35] 
have been developed to efficiently encode the full dynam¬ 
ics of DE/MG models featuring a single new degree of 
freedom. Linear structure formation in a very general 
class of DE/MG models (Horndeski scalar-tensor gravity 
[36, 37] and beyond [35, 38, 39]) can be fully described 
in this way using a small set of functions of time only. 
We can thus use this formulation to compare QS and full 
solutions in a general way in order to understand the ap¬ 
proximations properly. In what follows, we will use the 
term “dark energy” to refer to the models described by 
the EFT, i.e. to also most models of modified gravity. 

In this paper, we show that the domain of validity 
of the QS approximation is at best determined by the 
sound horizon of dark energy: even when the parameters 
of the DE models evolve slowly, the QS approximation 
can only be trusted inside the sound horizon. The scale 
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of this sound horizon can be calculated from an algebraic 
relationship between the parameters of the EFT or, al¬ 
ternatively, is a functional of the action for the DE model. 
We derive a condition determining when the QS approx¬ 
imation can be used safely to analyse a particular dataset 
and apply this result to parameters of current and future 
surveys. 


II. THE QUASI-STATIC SOLUTION 


We model the late universe as a mixture of pressure¬ 
less dust with density fraction Q m , representing both 
dark matter and baryons between which we do not dif¬ 
ferentiate, and the dark energy with density fraction 
Ode = 1 — Q m , which we assume comprises a single de¬ 
gree of freedom. The background expansion history is 
constrained by observations (e.g. ref. [40]) to be that cor¬ 
responding to a source with approximately constant pres¬ 
sure and no spatial curvature. In the linear perturbed 
Einstein equations describing the evolution of structure 
on largest scales, an extra scalar degree of freedom can 
be integrated out by solving the constraints. This pro¬ 
cedure, described in detail in ref. [34], reduces the exact 
perturbation equations into a second-order differential 
equation for the gravitational potential <f> coupled to the 
pressureless matter. The full form for such an equation 
for general DE models which do not contain operators vi¬ 
olating Lorentz invariance is given in [34, 41]. 2 However, 
for the purpose of this work, it is enough to write as a 
proxy an approximation valid at smaller scales, where it 
can be written as 
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where we use primes to denote derivatives w.r.t. e-folding 
time, S and v are the matter density contrast and velo¬ 
city potential in Newtonian gauge, H is the Hubble para¬ 
meter, a the scale factor, kn = k/aH is the comoving 
mode number in units of the cosmological horizon, while 
7 is a model-dependent friction coefficient. The sound 
speed c s > 0 is the speed at which high-frequency longit¬ 
udinal (scalar) perturbations in the DE propagate, while 


2 The beyond-Horndeski models of ref. [35] contain higher- 
derivative operators which cancel once the constraints are solved 
and therefore their dynamics reduce to (1) in the appropriate 
limit. Including true Lorentz-violating operators, brings higher- 
order spatial derivatives, fc 4 -terms, to eq. (1). These sort of 
terms would lead to a violation the r~ 2 law in gravity, which 
constrains them to lengthscales l < 1 //rri [42]. There must there¬ 
fore be a large hierarchy between the sound-speed and the A; 4 
terms. In such a case, including such suppressed fc 4 terms does 
not change any of the conclusions of this paper. 


M is some effective mass of the extra degree of freedom 
related to parameters of the model and the background 
expansion history (see e.g. [34]). The potential is coupled 
to matter through an effective Newton’s constant, p,s, 
and an equivalent coupling to the matter velocity poten¬ 
tial, fj, v . We take both of these to be constant, but time 
evolution does not add quantitative changes provided it 
be slow, with timescale H. The system is completed by 
the standard EMT conservation equations for dust 

S' = k 2 H Hv + 3$', Hv' = -$ , (2) 

where we have assumed that there is no gravitational slip, 
$ = VJ/. 3 The two gravitational potentials are always re¬ 
lated through a constraint rather than a dynamical equa¬ 
tion, and therefore generalising this result to models with 
non-vanishing slip does not change qualitatively the res¬ 
ults presented below: if the quasi-static approximation is 
good for <t>, it is also good for T. 

The coupling to the velocity /i v always appears when 
the Einstein (00) and (Oi) equations are used to eliminate 
the DE scalar perturbation and it does not in general dis¬ 
appear when the comoving-gauge density perturbation is 
used. Secondly, we stress that since the modification of 
gravity contains an additional scalar on top of the stand¬ 
ard scalar density perturbations in matter, the time de¬ 
rivatives of the potentials $ depend on two time-scales: 
(i) that of the pressureless collapse of the dust, which is 
just H and does not produce any oscillations since there 
is no pressure (ii) of the modihed-gravity scalar degree of 
freedom, which is related to its mass and the sound hori¬ 
zon and thus to the pressure support it provides. This is 
somewhat different than in ref. [44] , where the time-scale 
intrinsic to the scalar was not considered to determine 
the oscillations of the gravitational potentials. 

In principle, to solve the system (1-2), we should find 
the eigenmodes and their initial conditions. This is in 
general impossible. We instead can define the QS solu¬ 
tion by splitting the Newtonian potential into an oscil¬ 
lating part ip and the quasi-static part 4 >qs 

$ = ip + $ Q s , (3) 

choosing 4 >qs such that the oscillating part <p has no 
source in its equation of motion, £[<p\ = 0, with £ de¬ 
noting the homogeneous part of eq. (1). This gives us 
one decoupled equation for ip and the dependent eqs. (1- 
2). Since the equation for ip is source free, provided that 
its coefficients are appropriate, the oscillations ip can be 
expected to decay away over time , leaving us with only 
what we have called the QS solution. This is close but not 
exactly the same definition as is used when time deriv¬ 
atives are just neglected. We include extra information 


3 In the EFT, this sets both the parameters aj = = ®H = 0 

[43]. 
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which allows us to make a definite statement about the 
region of validity, inside of which the two solutions are 
the same. 

As we will see, there are some corrections to £ [p] which 
are suppressed by (c s /c#)~ 2 . These corrections will there¬ 
fore place an ultimate limit to the scales on which the QS 
limit can be considered to be valid. 


III. MAXIMUM DOMAIN OF VALIDITY 


In general, we cannot find an exact solution for the 
particular integral of (1). We can, however, obtain the 
Pade approximant to this solution, valid at least on some 
scales. Therefore, for the QS solution, we will be seeking 
an approximate solution valid inside some set of small 
scales to some accuracy. 

We start off by considering the typical evolution of 
cosmological-perturbation variables inside the cosmolo¬ 
gical horizon, ku ^ 1. In the usual case, we have 
Hv ~ 0(k]fS) by virtue of the continuity equation (2). 
This gives a natural hierarchy for these variables which 
will allow us to carry out an expansion. In what follows, 
we will therefore consider 

5 ~ O(e 0 ), Hv ~ 0(e), (4) 

identifying the order parameter e roughly with kff'. This 
also gives us $ ~ 0(e). By construction, our approx¬ 
imation always fails at the cosmological horizon, if not 
before. 

Given this hierarchy, we can choose an ansatz for the 
QS solution, which includes the contributions lowest- 
order in e and is compatible with eq. (1), 

= ~l(Ai + ^Sl m 6-^n m Hv. (5) 

The time-dependent coefficients A\,A 2 Bi are chosen in 
such a way that eq. (1) can be written as £\p\ =0 up 
to the relevant precision defined by e when <f> is replaced 
according to eq. (3). 

In the extreme quasi-static limit, kn -A oo, one 
typically neglects the A 2 and B\ terms, since they are 
suppressed by e. However, if we wish to include, for 
example, the corrections resulting from the mass M 2 , 
one needs to include terms such as A 2 and therefore also 
B\ since both of these corrections appear at the same 
order. Moreover, as we will see, including B\ allows 
us to estimate the domain of validity for the QS solution. 


as 

( 6 ) 
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with the J~, functions of time. According to how we have 
defined the QS solution, we need to pick such coefficients 
A i t 2 and B\ that all the 8F t = 0. The particular choices 
we require are: 

Ai = ps , (7) 

C S A 2 = -T^rn/U “ M 1 _ #" + #2 ~ ^ ) > 

CgBi = Cg fji V + n&(A - 7). 

This choice fixes all the freedom in ansatz (5) and 
removes the external source in the evolution equation for 
p, eq. ( 6 ), decoupling it from the matter perturbations. 
This choice also matches the standard choice for the 
leading-order terms of the QS solution. In this new 
description, p evolves with no source and, given certain 
additional conditions, p decays away and the QS solu¬ 
tion is eventually reached. We discuss these additional 
conditions in section IV. 

Let us now turn to the evolution equation for S, or 
equivalently, the equation for the growth factor, which 
is usually the equation of interest when the QS approx¬ 
imation is used. Using the QS ansatz given by (5) with 
the coefficient values given in eqs (7), we eliminate the 
matter velocity potential v from the matter conservation 
equations ( 2 ), obtaining 

S + ^2 + + eJ- 4 ^ 8 — (1 + eT, 5 ) 8 = S'[ 92 ] ( 8 ) 

where the source is determined by the deviation away 
from the QS solution, 

S[tp] = - (^kjj - < ^eiu, s fl m S J p + 3e ^2 + tp' + Sep”, 
and 


H 1 

-H'V 


c sk 2 H Fi = n m Q W (7 - 4)) + 0(c 2 ), 

c 2 s k 2 H T .5 = + 1 - ^ - 7 + 0 (c 2 ). 


One now substitutes the QS ansatz (3) into eqs (1), (2) 
and (3), keeping all terms involving p and replacing any 
time derivatives of 8 and v using the conservation equa¬ 
tions. The evolution equation (1) can then be rewritten 


In the small scale limit of eq. ( 8 ), kn —>• 00 , we re¬ 
cover the standard equation for the growth factor with 
sources which all involve the oscillating part of the poten¬ 
tial, p. Thus if indeed we can show that p decays away 
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at these scales, the QS solution for the growth factor will 
be valid once that has occurred. One should note that 
as the mode evolves from superhorizon to subhorizon the 
deviation from QS <p does contribute to the dynamics. 
Thus even if deep inside the sound horizon the QS rela¬ 
tionship between <f> and 8 is correct, it is not necessarily 
true that the standard relationship between the super¬ 
horizon curvature perturbation and the subhorizon 8 is 
maintained (see section IV for a detailed discussion). 

As we move away from the extremely small scales, 
corrections to the evolution appear. In particular, if 
M 2 H 2 , then the corrections due to this mass are as 
expected from the form of the coefficient of $ in eq. ( 1 ) 
and as is currently taken in f(R) gravity models [45]. 
There are also some additional corrections proportional 
to psH 2 which can only become relevant when the QS 
approximation in breaking down close to the cosmolo¬ 
gical horizon. Thus the naive QS estimate of the effect¬ 
ive Newton’s constant is essentially valid at cosmological 
subhorizon scales. 

Crucially, there are also corrections to the friction 
term, S'. The corrections are of order ps and are only 
suppressed by the sound horizon for the dark energy, 
c s kn rather than the cosmological horizon kn- Thus, 
at scales near the sound horizon, the response of the 
growth function to some particular value of the effective 
Newton’s constant is different than what would be 
expected from the naive guess, even if the quasi-static 
solution for <f> were valid. The typical sign of this 
correction is negative, thus the growth rate increases as 
the sound horizon is approached. This suppression by 
only the sound horizon is driven by the fact that the 
solution for the coefficient B\ contains terms of order 
c ” 1 and thus such corrections are missed when the dust 
velocity potential v is ignored. 

We should also add that similar corrections appear to 
the evolution equation for ip, eq. ( 6 ) at order e. The 
coefficient of the friction term is also modified by a term 
J : i/c 2 k'jj, thus leading to a similar decrease in the friction 
and therefore a relative enhancement of the amplitude of 
the oscillation ip. Thus the sound horizon also provides a 
limiting scale for the quasi-static approximation to hold 
for the potential $. 

We also remind the reader that in the derivation above 
we have neglected all time derivatives of the DE model 
parameters. If the time scale of the evolution is of order 
H , then including these terms does not qualitatively af¬ 
fect the results. However, if the model features a rapid 
time evolution of its parameters, e.g. with timescale H/e, 
then these corrections can be large and would dramatic¬ 
ally change the conclusions. Unsurprisingly, if the model 
features a rapid time evolution of its properties, one can¬ 
not hope that a quasi-static approximation can hold. 


IV. EVOLUTION OF OSCILLATIONS 

In the previous section, we demonstrated that, inside 
the dark-energy sound horizon, the quasi-static solution 
can in principle model the evolution of the full dark- 
energy-dust system. Choosing a scale well inside the 
sound horizon, c s kn 1 , is a necessary condition for 
the existence of this QS solution. It is not sufficient, 
however. In addition, the deviation from the QS solu¬ 
tion ip must decay faster than the QS solution, such that 
the QS solution is asymptotically reached after sound- 
horizon crossing. In this section, we will demonstrate 
the conditions required for this decay to occur. 

The definition of the QS solution we have proposed in 
section II was constructed in such a way that the variable 
ip — describing the deviation of the gravitational poten¬ 
tial from the QS solution — has no source in its evolution 
equation up to some precision. This decouples the evolu¬ 
tion of ip from the dust part of the system and therefore 
its behaviour can be considered in isolation. Taking the 
small-scale limit, kn —> oo, where the QS solution should 
be valid, if it is valid anywhere at all, the variable ip obeys 
the homogeneous part of eq. ( 6 ), 

tp + - l)Htp + c 2 k 2 ip = 0, (9) 

where we have re-expressed eq. ( 6 ) in conformal time r], 
with the overdot signifying a derivative w.r.t. it, while 
7i = aH is the conformal Hubble parameter. 

For any particular model with its appropriate back¬ 
ground, eq. (9) can be solved numerically. Here, we will 
only perform an estimate in a situation where closed- 
form solutions are available: assuming power-law evolu¬ 
tion where necessary. The conformal Hubble parameter 
evolves as H = 2/ |1 + 3wtot| V where w to t ^ —1/3 is 
the total equation of state for the universe. We are 
also going to assume that the dark-energy sound speed 
evolves as a power law with redshift, c 2 = c^ 0 (l + z)~ p = 
c sO r ? 2p ^ 1+3lUt °^ • We note here that for p < —1, the co¬ 
moving sound horizon would actually be shrinking dur¬ 
ing matter domination and therefore the modes would be 
leaving the sound horizon. For any particular mode, the 
super-sound-horizon corrections to eq. (9) would eventu¬ 
ally take over and the QS limit would be violated com¬ 
pletely by the results of section III. 

The friction term 7 for general scalar-tensor models is 
determined by the action and given in full in ref. [34]. 
In MG models, where the kineticity is small compared 
to the braiding, qk <C Ob (e.g. f(R)), 7 ~ 3; on 

the other hand, in perfect-fluid DE models, such as 
k-essence/quintessence 7 « 4. This value is model- 
dependent, so we will only assume it is a constant during 
the era being considered. We can then rewrite eq. (9) in 
(one of the) standard forms of the Bessel equation, the 
solutions of which are oscillatory with a decaying envel¬ 
ope. Expressing the result as a function of redshift, the 
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envelope evolves as ip env oc (1 + z) v with decay exponent 



7~1 

2 


( 10 ) 


when inside the sound horizon. Note that this is inde¬ 
pendent of Wtot and is also a good approximation for the 
ACDM background. So if we take 7 = 3, then p > —4 
for a decaying ip. Thus the sound speed must decrease 
rapidly for p to grow. Beware that the comoving cos¬ 
mological horizon shrinks during acceleration (z < 0.65 
in ACDM with fl m = 0.31) and the modes begin to exit 
the sound horizon, potentially invalidating the QS limit 
in any case. 

In eq. (9) we have neglected the mass term M 2 a 2 . 
For example in f(R) models, this mass can be larger 
than the sound-speed term even on scales inside the 
sound horizon. On these scales, the evolution of <p env 
is determined by the evolution of the mass, but the 
solution goes through similarly provided it is a power 
law in (1 + z). The exponent v will change as the mode 
crosses the scale k(l + z) ~ M/c s . 


The conditions we have derived refer to the QS rela¬ 
tionship between the gravitational potential and the mat¬ 
ter density perturbation, as given by eq. (5). A separate 
issue is the relationship between the initial superhorizon 
curvature perturbation, as set up by inflation, and the 
amplitude of the mode inside the horizon. Every mode 
starts off superhorizon and will have evolved through 
scales where the QS limit fails, even if the scales probed 
by a the survey are in the deep QS regime. The valid¬ 
ity of the QS solution guarantees the evolution for the 
growth factor from one redshift to another, but does not 
necessarily ensure that the amplitude will not have been 
modified previously by the interaction between the DE 
degree of freedom and the matter on superhorizon scales 
and while the deviation from QS is dying off. Eq. (8) 
contains a source involving ip for a large fraction of the 
evolution history of all modes. This will end up having 
an impact on the shape and normalisation of the mat¬ 
ter power spectrum inside the sound horizon which in a 
generically scale-dependent manner. See ref. [46] for an 
explicit example. The only way to suppress this effect is 
to ensure that for the modes being observed, the modi¬ 
fications of gravity were insignificant at the time that the 
modes were outside of their sound horizon. 

We discuss the impact of the above results on current 
and upcoming surveys in the next section. 


V. APPLICATION TO SURVEYS 

In the preceding, we have built up a picture where, 
outside the dark-energy sound horizon, the configuration 
of the dark energy degree of freedom is significantly dif¬ 
ferent from the QS solution. Provided that the comoving 


sound horizon be growing, the mode will cross it. The 
deviation from the QS solution will then decay away with 
the exponent given by eq. (10), eventually becoming neg¬ 
ligible given some required theoretical precision. We will 
ignore the possible temporary effect of this decaying de¬ 
viation on the density perturbation in this section. 

Given the above, we can estimate the minimum dark- 
energy sound speed which is necessary for the QS approx¬ 
imation to be appropriate in the analysis of observations 
given particular survey parameters. We will make the 
assumption that the QS limit of the dark-energy model 
is very simple, with no rapid time evolution of its prop¬ 
erties or scale dependence. This is typical of the models 
which are considered in data analysis (e.g. ref. [47]). We 
assume that 


• the cosmological background is exactly ACDM with 
parameters approximating the Planck 2015 results, 
Sl m o = 0.31 and z eq = 3371 [40]; 

• the sound speed of dark energy c s is constant, 
p = 0; we take a model with 7 = 3, i.e. with the 
exponent of the decay of oscillations of eq. (10) con¬ 
stant, v = 1; 

• the value of the effective Newton’s constant in the 
QS solution is in principle redshift dependent and 
given by /r 5 (z); 4 

• the survey uses some longest mode fc m i n = 27r/r max 
for their DE-related analysis using a sample of 
galaxies within a redshift bin centred on z surv . 


We then require that, for the longest mode fc m i n , the 
oscillating part of the potential constitutes no more than 
fraction e of the QS value, i.e. the DE sound speed needs 
to be large enough to satisfy 


<p(z *) 


^QS (Zsurv) 


1 


1 + 


5, £ 


( 11 ) 


where z* is the redshift at which the mode k m - ln crossed 
the DE sound horizon, obtained by solving c s k m i n = 
Ti(z*). The deviation from the QS solution at sound- 
horizon crossing can be roughly estimated by assuming 
that above its sound horizon, the dark energy will behave 
roughly as dust (since there is no pressure support), so 
one would expect (Q m + IIde)^- On the other 

hand, the QS result is given by kjj^Qs = 

Thus 


< p ( z *) = - 4>qs 


~Vs 1 (v5~ 1 - $QS(^*) • 

( 12 ) 


4 Allowing for redshift dependence of fis does not change A± in 
eq. (7), which is all that is necessary for the discussion here. 
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Let us now parameterize the effective Newton’s constant 
as [ 20 ] 


M<5 — 1 + no(l + 2 ) s (13) 


and assume that fis is not too different from one, so that 
the potential is approximately constant during matter 
domination. Also, we assume that the dark energy is not 
of tracking type and therefore £Ibe{z*) D m (z*). In 
order to obtain an analytic solution, we will also approx¬ 
imate the Hubble parameter as H « Hoy/£l m o(l + z), 
which is a reasonable approximation during matter dom¬ 
ination. Putting the above together, we obtain from con¬ 
dition ( 11 ) 


> Vn 


mOl 


Ho ( Mo \ 


e / 


(i 


v/{a+v) 


(14) 


This should be thought of as an order-of-magnitude ap¬ 
proximation, which should be compared with numerical 
results in table I. We can easily see that if the QS solution 
is to be good enough, we need higher sound speeds for 
high-redshift surveys, large maximum lengthscales sur¬ 
veyed, larger deviations from /ro = 0 and more slowly 
evolving effective Newton’s constants, just as one would 
naively expect. The minimum sound speed is weakly de¬ 
pendent on the survey redshift, but strongly on the low¬ 
est mode fc m i n being used in the analysis. Eq. (14) is the 
main result of this paper, allowing for a simple estimate 
as to whether the QS approximation can be used for a 
given dataset and DE model. We remind the reader that 
this approximation is valid for simple DE models with 
constant sound speed and a monotonic evolution of the 
effective Newton’s constant. 

In table I, we have summarised a selection of surveys 
and the constraints on the sound speed given the require¬ 
ment of quasi-staticity. It shows that today’s galaxy sur¬ 
veys are well inside the region where the QS limit is valid 
for our simple models with sound speed c s > 0.02. There 
is an interesting effect whereby, for the surveys centred on 
lower redshifts, such as BOSS, the QS approximation is 
better than expected since the comoving horizon shrinks 
for 2 < 0.65 and therefore the modes fc min entered the 
sound horizon at higher redshift and are already leav¬ 
ing it at the time of the survey. Once Euclid data are 
available, the quasi-static approximation will only be ap¬ 
propriate for dark-energy models with sound speeds close 
to that of light. 

On the other hand, the CMB anisotropies already 
probe very large scales, since the Hubble horizon at 
recombination is approximately one degree on the sky 
today. This means that if the DE was playing any role 
at all at recombination (e.g. in Early Dark Energy mod¬ 
els [48, 49]) the quasi-static limit cannot be consistently 
used there, even for models with c s = 1. Whether it is 
consistent to use the QS approximation for late time DE 
in CMB data is, however, subtle. One must ask both 


which range of scales k contribute to the kernel of Cf. for 
the particular multipole of interest £ and whether at the 
redshift at which they contribute, the DE perturbations 
are at all important. Only if they are, then one must con¬ 
sider whether the perturbations of the DE and whether 
they are already quasistatic at those scales. 

The DE perturbations have two main effects on the 
CMB: a modification of the integrated Sachs-Wolf effect 
(ISW) and the lensing of the CMB. the full the ISW effect 
has its kernel peaked at approximately k = 0.004 Mpc -1 
for £ = 20 and k = 0.0005 Mpc -1 for £ = 2 and it is sens¬ 
itive to exactly the late times. Thus one should not use 
the QS solution at all (see e.g. [50, fig. 4]). Nonetheless, 
the impact on parameter estimation is limited because of 
cosmic variance at such scales. 

On the other hand, Planck is sensitive to CMB lensing 
in multipoles 40 < £ < 400 [51]. Using ACDM results as 
a proxy, for £ = 40, half of the signal is from z < 1; even 
at £ = 400, half is still from z < 2 [52], These are the 
times at which DE dominates and therefore here also one 
must be certain that the DE perturbations are calculated 
properly at the scales of concern. For £ = 40, the kernel 
is peaked at k = 0.01 Mpc -1 , while for £ = 400, the 
peak is somewhat below k = 0.1 Mpc -1 . This implies 
that one should not use the QS approximation in CMB 
lensing predictions whenever c s < 0.1. 


VI. DISCUSSION AND CONCLUSIONS 

Adopting the QS approximation essentially means that 
the DE degree of freedom with a finite propagation speed 
is replaced by a modification of the constraints of gen¬ 
eral relativity, i.e. with instantaneous propagation. The 
advantage is that the dynamics are much simpler, insens¬ 
itive to DE initial conditions and expressible in terms of 
simple parameters. However, this approximation must 
fail whenever causality places a limit, i.e. beyond the dark 
energy sound horizon. We have shown that for current 
galaxy surveys, the QS limit is sufficient in the analysis 
whenever the DE model has sound speed c s > 0.02. As 
we enter the era of much wider surveys, such as Euclid, 
this approximation becomes inappropriate even for mod¬ 
els with sound speeds close to those of light and could 
result in misleading constraints, both because the sound 
horizon is too close to the scales involved in the analysis 
and because the closeness of the cosmological horizon in¬ 
validates the QS approximation in any case. Moreover, 
the CMB probes the universe at the largest scales already 
today. Our result implies that the QS limit is never valid 
for early dark energy, good enough for CMB lensing if 
c s > 0.1 and always inappropriate for ISW since scales 
close to the cosmological horizon are relevant. 

On the other hand, most of the well-studied models 
of dark energy, e.g quintessence and f(R ), have sound 
speed equal to that of light. The QS predictions for ISW 
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Minimal sound speed for good QS limit 


Survey: 

BOSS 

CFHTLenS 

DES 

Euclid 

•2-surv 

0.57 

0.75 

0.6 

2 

fomin/h MpC 

0.04 

0.2 

0.01 

0.007 

&min / ~L ( ^surv ) 

140 

700 

50 

20 

Cs,min@(£/ M0 — 0.01, S 

= 1) 0.02 

0.003 

0.04 

0.1 

Cs,min@(£/ M0 — 0.01, S 

= 3) 0.009 

0.002 

0.02 

0.06 


Table I. Minimal dark energy sound speed c a which allows for 
the use of the quasi-static approximation for given example 
survey parameters assuming Cs be constant. Current surveys 
can be safely analysed using the QS solution for c s > 0.02 
with accuracy better than e/po = 1%- Future surveys will 
probe much closer to the cosmological horizon and this sort 
of accuracy can only be achieved for Cs > 0.1, depending 
on the largest modes included in the analysis. When modes 
comparable to the horizon are included, the QS limit should 
not be used at all. 

BOSS assumes largest scales used in redshift-space-distortion 
measurements of r max = 152/i _1 Mpc [47]. For CFHTLenS, 
we assume largest angular scales available of 0 max = 100' 
[53], while for DES — that scales 3 times smaller than the 
maximum observed in the 5000 sq. degree survey will be used, 
i.e. 23° [54]. Maximum Euclid scale with significant data 
chosen as in ref. [55], i.e. 40°, which is approximately 1/3 of 
the maximum angular scale available. Actual precision of QS 
for Euclid is worse since cosmological horizon corrections are 
already significant at those scales.We employ a full numerical 
computation rather than the approximate result (14). 

in these models are in all likelihood somewhat off the 
full solution, but we readily admit that cosmic variance 
at those scales means that the biasing of parameter es¬ 
timation is small. For other probes, using the QS ap¬ 
proximation in these particular models is good enough 
for current surveys, as shown in the literature previously. 
Thus the effects such as a scale-dependent Newton’s con¬ 
stant in f(R) gravity when M 2 ^ k 2 /a 2 do exist and 
might possibly be observed [45]. 

We have provided condition (14) as a simple test allow¬ 
ing us to determine whether the use of the QS approxim¬ 
ation could at all be valid when we use observational data 
to test models of gravity beyond the simplest ones such 
as f(R) and quintessence — for example, those defined 
by their EFT parameters. This test should be treated as 
a best-case scenario, which would be violated by a rapid 
evolution of any of the parameters. One should also bear 
in mind that there may appear a scale-dependent modi¬ 
fication to the amplitude of the matter power spectrum 
resulting from some dynamics during DE sound-horizon 
crossing, even for modes which are well inside their QS re¬ 
gime during observations (see ref. [46] ). As we map larger 
scales in the near future with Euclid, LSST and the SKA, 
an EFT-like formulation including all the time-domain 
behaviour will increasingly become necessary to connect 


consistently initial conditions at large scales (e.g. as given 
by the the cosmic microwave background) with late-time 
observables. Fully implementing the dynamics need not 
necessarily lead to a degradation of constraints, but has 
been demonstrated to sometimes improve them [56]. 
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